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(and the opposite in sign excess charge at the surface) 
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caused by the influence of gravitational forces is determined. 
in 

Up to now an old problem of connecting rotatoty and magnetic charac- 

> ; teristic heavenly bodies has been of gread interest. One of the aspects of 

(S| ! the above problem is a question concerning quasi-neutrality disturbance of 



the substance of the body, (neutrality which on average in scale is greater 

than the lattice parameter or Debay radius because of gravitational forces. 

In addition to condensing the substance of the body these forces also create 

^ an excess of a positive charge in the center of a heavenly body and an excess 
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of a negative charge on periphery of a body, (gravitational forces affect more 
strongly heavy nuclei than light electrons). The electric field caused by this 
effect results in the compression of the electric component of the substance 
of a heavenly body. 

Rotating around its own axis that body owing to redistribution of its 
charge must obtain magnetic moment. Thus far the scientists have tried 
to explain (by the mechanism under discussion) the magnetic properties of 
heavenly bodies. The aim of this paper is to prove hopelessness of such 
attempts due to a negligible value of the effect of charge redistribution which 
corresponds to parameter a: 






a = —^ « 1(T 36 , (1) 
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Here G - a gravitation constant, m p and e - the mass and charge of a proton. 
1. Let's start with elementary evaluations directly concerning a solid (crys- 
talline) state of a substance. In this case redistribution of a charge is the 
displacement of a nucleus relative to the center of "Vigner-Zeitc cell" and 
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leads to the appearance of dipole moment and polarization [1]. It is obvious 
that both the electric and gravitational forces which act on the displacement 
of a nucleus must be in balance. Hence 

ZeE = -GM 2 v|dx / n p (x / )/|x-x / |, (2) 

Here Ze and M - a charge and mass of a nucleus, E - the vector of an 
electric field, n p - the concentration of nuclei. Taking into account equation 
div~E = 4tt5p (5p - an excess of a charge), one may find: 

S P = (^f a Pp ~ a P P - ( 3 ) 

Here A - the mass of a nucleus divided by the mass of a proton, p p - the 
density of an electric charge of nuclei. 

Simple dimension considerations are in favour of evaluation (3). 

Dimensionless relation bpj p p must be proportional to G (or more exactly 
to a) which follows from the application of perturbation theory on the grav- 
itation interaction. Corresponding coefficient of proportionality may depend 
on dimensionless parameters of substance - z, the ratio of masses of a proton 
and an electron; the ratio of Coulomb energy of a single particle to the largest 
of the following values: Ep and T [Ep - Ferme energy and T - temperature 
). Now you can see for yourself that for astronomic bodies the above pa- 
rameters differ from 1 not more than 100 - 1000 times. Therefore they can't 
considerably influence eq. (3), (because 5p/p p , depends on G linearly). The 
negative charge, which compensate (3) and is localized on the surface of a 
body, equals 

Q = — dx.5p ~ aQ p , (3') 

Here Q p - the total charge of the nuclei of a body. The polarization P = 

Zen p 5 (5 - the shift of a nucleus) is equal to — E/47T, because induction 

D = owing to the absence of external charges. As a result the elementary 

calculation using eq. (3) gives the ratio of the magnetic moment of a body 

to its mechanical moment. 

~ -a . (4) 

m p c 

Owing to the infininitesimal of parameter a (see (1)) the above mentioned 
values are very small and the effect of charge redistribution under discussion 
can't have direct observing demonstrations. 



Suffice it to say, that for the Earth (the mass ~ 10 27 g., the radius ~ 10 9 sm) 
the magnitude of the surface charge (3') corresponds to one electron per 1 
m 2 of the surface. 

2. A more rigorous derivation of eq. (2) for a solid state of a substance 
is based on the selection of the part system from the total energy depending 
on the shift of nuclei 5k (k - the number of a nucleus) and minimization of 
this part over 5k accompanied by a change of: 5k — ► Pk/(Zen p ). Besides 
the sums for the grate may be replaced by integrals: 

]T — ► f d-Kn p , p k — ► P(x). 
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Let's start with gravitational energy for nuclear interaction : 

&gr = ~ — Z^ l X k — x k'| 

1 k,k' 

Substituting Xk — > Xk + 5^ and factorizing over 5 up to the first order in- 
clusive, one finds: 

GM 2 f , rn , x _ r (hdnJ-x!) , 
2^ = -^/dx{P(x)v/^f^i}. (5) 

Coulomb energy of interaction of nuclei and electrons can be written as: 

E|xk-x / k |- 1 -Ze 2 /dxn(x)E 1 ^ "' ' ' 
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^ c = ^^E|xk-x / k |- 1 -Ze 2 / G ?xn(x)E|xk-x|- 1 +- o 
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here n ~ Zn p - the concentration of electrons. In repeating the same consid- 
erations which resulted in (5) we must have in mind two circumstances: 

1. Now we must factorize 5 (or P) to the second degree inclusive, (as a is 
very small, see (1)). 

2. Factorization with respect to 5 is impossible in case of interaction with 
electrons of the same cell because of the necessity to take into account the 
contribution of region r < 5 and this case must be considered exactly. 

It is represented in the second term of E C1 in which one must select an 
integral over the cell volume: 

dxn(x) 
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|x-(5|' 

multiplied by the total number of nuclei. It leads to the following expression 
for the part depending on P: 



In the remaining part of E c after factorizing over S an ordinary term of dipole- 
dipole interaction arises, the interaction of different cells being reduced to it: 

E m = f /^ ( WW) _ 3(p(x),p(x-x'))(pw,(x-x0) 

J J V 2|x-x'| 3 2|x-x'| 5 ' 

(Terms of E c which are linear with respect to P disappear owing to the neu- 
trality of a system). The equation for EJ?> may be written as (see Enclosure): 

£< 2 > = !/dx[P 2 (x)-3(P(x),V)^»] 

For that reason the sum of E^ and E^ can be written as: 

£ c = 27r|dxP 1 2 , (6) 



here Pi = ^^p . the longitudinal part of vector P (it is in fact in (5)). 
Minimization of sums (5) and (6) over Pi (if E = —4ttP) brings us back to 
eq. (2). 

Let's emphasize that in the last equation there is no correction for the 
distinction of the functioning field from the average one [1] . That correction 
would appear if we placed the nucleus considered - for which the balance of 
forces is written - into empty space (2). In fact we have the interaction of 
a nucleus with electrons of its cell (cavity). It is described by E c and the 
required contribution — 47rP/3 to intensity because |jy = — E. To complete 
the proof of eq. (2) one should make sure that the parts of the system energy 
not considered above (to be more exact their parts depending on P) do not 
influence the result. First of all it is true for the electronic component of 
energy - kinetic, exchange, correlation [2]. 

In case of a strongly compressed substance - this case is more interesting 
for heavenly bodies having a solid state substance - kinetic energy of free 
electron gas is an important factor (other components of energy are very 
small compared with Coulomb energy taken into account above). 

Decomposition of this energy over the nuclear shift relative to the center 
of a cell: 

SE kin ~ [dxnSU, 5U = Ze\- - — |— ) 
J r |r + o\ 

leads to the zero result owing to the spherical symmetry n and a known 
decomposition of Coulomb term into a series over Lezhandr polynomial. As 
regards the energy of nuclei in a solid state at low temperatures one should 



consider only zero energy of nuclear oscillation which is equal to |ftu>o for 
one nucleus (wo - the oscillation frequency). Hence only part wq depending 
considerably on shift 5 or P could affect the proof given above. But in case 
of an exceptional small value of S (5 ~ aR - see point 1, R - the radius of a 
body) the shift of the oscillation center of a nucleus in the cells model does 
not influence the oscillation frequency at all. A change in energy in shifting 
a nucleus by 5 r is equal to l 2 e 2 5 2 /(2R: i ) (where R - the radius of a cell) and 
in substituting 5 r — ► 5 r + 5 (the shift of the oscillation center) the square 
of frequency as the second derivative of energy with respect to S r does not 
change at all. 

3. Eq. (3) is in fact of universal nature and it is valid irrespective of 
the state of a substance. We will utilize a method of functional density [3], 
writing down the free energy of a system (in general case temperature T is 
not equal to zero) as a function of the density of electrons n and nuclei n p . 

F{n, n p ] = F + E 3C + F xc - ji f dxn - p p J dxn p , (7) 

here the first term conforms to free electronic and nuclear gases, the second 
one - Coulomb and gravitation energy of their interaction in approximation 
of a self - consistent field, the third term corresponds to exchange and corre- 
lation effects. The latter two terms correspond to Lagranzh factors and show 
conservation of total numbers of electrons and nuclei. Minimum F for n and 
n p defines equilibrium distributions of these quantities. Writing down 

GM 2 r dxdx' . e 2 r dxdx! . w . .. 

E sc = — J |x _ x/| n p n p + - J |x _ x/| (n - zn p ){n - zn p ), 

we obtain conditions of minimum (7) for n and n p : 

5(F + F XC ) 2 

A — -— - — = 47re (n - zn p ) [7) 

dn(x) 

A S(Fo + F ) = _ 47rze 2 (n _ j _ 4nGM 2 ^ 

on p (x) 

If we could omit the left part of eq. (7"), then, taking into account that 
Sp = e(zn p — n), p p — zen p , we would obtain eq. (3), and if (7") is substi- 
tuted into {!') and if appropriate simplifications are made then (7') will fit 
Chandrasekar equations which define equilibrium configuration of electrons 
and nuclei. Let's emphasize that the left part of that equation is defined by 
the lightest particles - electrons, the right part - contains only gravitation 



quantities after the above substitution of (7") into (7') although gravitation 
does not act on the electrons directly. The electrons are affected by the 
electric field investigated in this paper, which only quantitatively (see (3)) 
coincides with a gravitational field. 

Thus eq. (7') and (7") are rewritten as: 

Sp = aPp(1 + a) -^ = MEl^., (8) 

where here and below F = Fq + F xc . 

The appearance of a either keeps the magnitude of Sp/ p p constant, or 
decreases this relation. The single case when it may increase considerably, 
represents exceptional nearness of a and — 1. But this is practically unbeliev- 
able. Furthermore, we may think that a << 1 . Let's illustrate it by using 
two examples [2]. For both examples it is supposed that the electron gas is 
strongly compressed and degenerated so that the corresponding contribution 
to F is ~ — fdxn 5 / 3 . In the first example nuclei are localized in nodes of 

TO J r 

a grate and the energy of their zero oscillations represent 8F ~ 4= / dxn 3 / 2 
(see above). Then for a in (8) we find: 



^Z(a nV3)-V2 « lt 

where m - the mass of the electron, ao = % /me 2 - it's "Bohr-radius". This 
smallness is connected with inequalites: m/M « 1 and aon 1 / 3 >> 1 in a 
compressed substance. 

The second example represents a weak non ideal nuclear " Boltzman-system' 
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for which F{n p ) ~ — e / dx^fe- is " Debay-Hukkel" - correlation correction. So 
for a one finds : -p7=^ — /tt: << 1 because the condition of weak non-perfection 

[yZaan*-' 6 ) 



is T » eV/ 3 . 

4. In conclusion let's go back to the question about the minimum of (7) 
in connection with disturbance of local electric neutrality of the system. 

Note that the above violation is typical of the crystalline state of a sub- 
stance in the absence of gravitation forces, which is evident and the electrons 
are not localized in contrast to nuclei which are localized at the point. 

It is important that this violation is not described by minimum of F, as 
a point in functional space in which its functional derivative disappears. In 
this case we deal with a boundary minimum reached when the parameter 



which defines the length of nuclear localization tends to its limiting value 
which equals zero. 

Let's consider the second clear Coulomb component of E sc which should 
be broken into neutral and spherical on the whole "Vigner-Zeitc" cells with 
a nucleus in the center. In such a model the radius of the cell is R and 
the nucleus is spread in sphere with radius p. This model represents energy: 
E sc = - jjj ^- 2X (t^+ x X\)^ wnere x = p/R- At x = this equation gives a 
well-known binding energy of the grate: 10 d • It is evidentl that maximum 
of \E SC \ is in fact reached on the boundary of permissible area (at x — and 
fixed R). 

Enclosure 
Initial expression for E^ can be written as: 

^ 2 ) = /dxP i (x)^ ? (k)P J (x), 

where k = — zV, V act on Pj(x), and 

i r d\ , „ _ XiXj ^ 



K,M = 1 -J^(6 i3 -3^)e 



This integral satisfies evident condition K^ = and hence it may be repre- 
sented as 

The expression for E)r mentioned in this paper gives the calculation that is 
not complicated but uwkward. 



In conclusion we would like to note critical debate on this questions with 
Vasilyev B.V., Grigoryev V.I. and Maximov V.I.. 
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